{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import math\n",
    "import matplotlib.pyplot as plt\n",
    "from scipy.fftpack import fft,ifft\n",
    "import pandas as pd\n",
    "from scipy.optimize import minimize"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "data081 = pd.read_csv(\"energy_2008_1.csv\")  \n",
    "x081 = data081['x']\n",
    "y081 = data081['y']\n",
    "xerr081 = data081['xerr']\n",
    "yerr081 = data081['yerr']\n",
    "\n",
    "data101 = pd.read_csv(\"energy_2010_1.csv\")  \n",
    "x101 = data101['x']\n",
    "y101 = data101['y']\n",
    "xerr101 = data101['xerr']\n",
    "yerr101 = data101['yerr']\n",
    "\n",
    "data102 = pd.read_csv(\"energy_2010_2.csv\")  \n",
    "x102 = data102['x']\n",
    "y102 = data102['y']\n",
    "xerr102 = data102['xerr']\n",
    "yerr102 = data102['yerr']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAn8AAAHvCAYAAAAlwLklAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3Xt41OWZ//H3DYQkIKc0wPZgxQpSD8EGUsUD7aLU4mlbsZBWWzR2EbvVsrVSty34U3StrXYXwS0Vdk11sS3EU6Va6DarRSy1JcHYohy0CwpVCARIlCTkcP/+mGRIwmQymflOJofP67rmIt/DfL830Kve3M9zP4+5OyIiIiLSN/RLdQAiIiIi0nWU/ImIiIj0IUr+RERERPoQJX8iIiIifYiSPxEREZE+RMmfiIiISB+i5E9ERESkD1HyJyIiItKHKPkTERER6UMGpDqA7iw7O9vHjBmT6jBEREREOlRSUrLf3Ud2dJ+SvyjGjBnDpk2bUh2GiIiISIfMbFcs92nYV0RERKQPUfInIiIi0oco+RMRERHpQ5T8iYiIiPQhSv5ERERE+hB1+4qIiEhUjY2N7N+/n0OHDtHQ0JDqcPqk/v37M3z4cLKzs+nXL7HanZI/ERERiWr37t2YGWPGjCEtLQ0zS3VIfYq7U1dXx969e9m9ezcf/ehHE3qehn1FREQkqvfff58Pf/jDDBw4UIlfCpgZAwcO5MMf/jDvv/9+ws9T8heBmV1hZssPHz6c6lBERES6hUSHGiVxQf0d6G8yAndf4+43DBs2LNWhiIiI9Br5D20k/6GNqQ6jz1PyJyIiIhKwdevWMWXKlJjv37JlCx//+Mepra1NYlQhSv5ERESk17jttts444wzGDp0KB/60IeYM2cOFRUVre559NFHOeWUUxg0aBDnnHMOJSUlra5v2rSJs88+m0GDBnHKKaewcuXKVtdffvllPvWpTzF8+HBGjx7NV77yFQ4cOBC+7u5885vf5M477wyfKy4u5qKLLiIrKwszY/fu3a2eecYZZ5Cbm8uDDz4Y1B9Fu5T8iYiISJc4Wt/Ia3+rZF9VTdLe0b9/f1auXMmBAwcoKytj9+7dFBQUhK9v2LCBr33tayxbtoyDBw9y1VVXcemll1JZWQnA4cOHueSSS7jqqqs4ePAgP/nJT7jxxhvZuDE0XN3Q0MDll1/O+eefT3l5Oa+//jp/+9vf+MY3vhF+x29+8xuOHj3K1KlTw+cGDx7M7Nmz+e///u92Y7/++utZunQpjY2NQf+xtKLkT0RERLrEnkPVVNXWs6T4jaS945577iE3N5e0tDRGjhzJTTfdxAsvvBC+vmLFCmbMmMHFF19Meno68+fPJz09naeeegqAJ598kszMTL797W+Tnp7OZz7zGa688kqWL18OhJLD/fv3U1BQQFpaGllZWcyaNYuysrLwO55++mmmTZvWqjN68uTJXHvttZxxxhntxv6pT32Kd999l1deeSXgP5XWtM6fiIiIBCpSU8fR+kb2VYXmsz328i627DnMwAHHalCr5p6blFiKi4uZMGFC+LisrIzrrrsufGxm5ObmhpO3srIyJk6c2CpxmzhxYrhil5WVxdy5c1mxYgV33303lZWV/OIXv+DKK68M319aWso111zT6VjT09MZN24cpaWlTJw4sdPfj5UqfyIiIpJ0ew5VHzvwNsdJ8sQTT7BixQoeeOCB8LmqqiraruYxfPjw8LBvR9cBZs6cyVNPPcXgwYMZNWoU/fr14zvf+U74+sGDBxk6dGhcMQ8dOvS4OYpBU+VPurXyI+XMXz+f+z99P9mZ2akOR0REYtC2irevsoYpP3w+fOxAZXUdS6/OZdSQjKTEUFRUxNy5c3nmmWdaVdGGDBlC23V8Dx06xCmnnBK+vnPnzuOuNydzO3bs4JJLLmH58uVcc801HDlyhPnz5zN9+nTWr18PwIgRI1oli51RWVlJVlZWXN+NlSp/0q395NWfULq3lGVly1IdioiIxGlJ8Q4a3Vuda3BP2ty/wsJC5s6dy5o1a1o1XQCcddZZlJaWho/dnVdeeYWzzjorfH3z5s2tvrN58+bw9bKyMkaMGMF1111HWloaw4YN4+abb+bFF18MJ5W5ubm89tprnY776NGj7Nixg9zc3E5/tzOU/Em3UrC2gIK1BUx4ZAI5j+SwettqHGf1ttXkPJLDhEcmULC2oOMHiYhIt1H61iHqGlonf3UNTumug4G/a8mSJdx6662sW7eO888//7jrc+bM4cknn6S4uJijR4/yox/9iJqamvCcvSuvvJIjR45w3333cfToUYqLi3nyySe54YYbAMjLy+Pw4cOsXLmShoYGqqqqePDBB/nYxz4WHi7+/Oc/T3Fxcav3NjY2UlNTE17Hr7a2lpqamladvevXr2f06NFK/qRvysnOISv9WNnbMLIyspiQPSHKt0REpDt6bt4Udt57GeecnMU5J2ex897L2HnvZTw3L/ZFkGM1b948KisrmTp1KieccEL40+yCCy7gxz/+MXPmzGHYsGGsXr2a5557LjysO3z4cJ577jmKiooYNmwYc+bM4Sc/+Qnnnhsayh4zZgxPPPEEDzzwAB/4wAcYM2YMu3bt4pe//GX4HZ/97GcZMGBAqy7j9evXk5mZycc//nEAxo4dS2ZmZnioGODhhx/m5ptvTvpWeuZtyrByTF5enm/atCnVYfRZizYuomh7EUao42rm+JksnLwwxVGJiPQ9r7/+OqeddlrCz2nuAk5WZ293snbtWu65555WyV00r732GldeeSWvvvoq6enp7d4X7e/CzErcPa+jd6nhQ7qtipoK8sfnM/PUmRRtL2J/9f5UhyQiIhKT6dOnM3369JjvP/3009m2bVsSIzpGyZ+kRCxdvIunLg7/vGDygq4KTUREkqQvVPx6As35k5RQF6+IiEhqqPInSdeyO7dkbwnOsXmmq7etZvW21RjGpNGTKJxeGNc7tB6giIhIbFT5ky6VrC5eVRJFRERio8qfJF3bat6ijYt4fPvjpPVPo66hjmknTet0F29zNTFZlcT2qMIoIiI9nSp/0uUqaiqYNX4WP7v0Z8waP4sD1QfiflZXrweoCqOIiPR0WucvCq3z1zO0rSQmsh5gpMpewdqC4yqMzQzj1Wtfjfp9EZGeLqh1/ii8LPRrwbOJP6uPCmKdP1X+pMfrTCWx/Eg51629rt01A9ur7MVaYcz/VT4le0siVgY7ereIiEhXUOUvClX+ep+7/nAXRduKWlUHY63stbfjSKzfv3D1hZRXlzNr/CztVCIiPYoqf523ZcsWZsyYwV/+8hfS0tI6vL+8vJy8vDxKSkrIzm5/5EiVvyQxsyvMbPnhw4dTHYokqGBtAQVrC5jwyARyHslh9bbVOM7qbavJeSSHCY+EKnexVPaadxwpuqLouApjtO83v7u8uhzguHd3RBVDEZHY3XbbbZxxxhkMHTqUD33oQ8yZM4eKiopW9zz66KOccsopDBo0iHPOOYeSkpLwterqambOnMm4cePo168fd99993HvOHLkCNdffz0jRoxg+PDhfPWrX6W6urrVPbfeeiu33XZbOPErKyvjkksu4YMf/CBmxoYNG1rdP3LkSK6++mruvPPOoP4o2qXkLwJ3X+PuNwwbNizVoUhAoiVnhdMLeeyyx7jopIvC1wCmnTSNlZetDH9n8dTFLJi8gPFZ41kweUF4B5KW3zeMgf0Htvp+ok0p7Q1FKykUkR7jvnFwxzDYtSH0uWNY6HPfuMBf1b9/f1auXMmBAwcoKytj9+7dFBQcW292w4YNfO1rX2PZsmUcPHiQq666iksvvZTKykoAzIzzzjuP5cuXc/bZZ0d8x7x589i6dStbt25l+/btvP7669xyyy3h69u2beOll17ii1/8YvjcwIEDmTFjBs8880y7sV9//fUUFhaGY0kWJX/SqxVOL+wwOWsWrbIXi/bmHsby7raJXLRq5aSVkwB1HotID/L+vs6dT8A999xDbm4uaWlpjBw5kptuuokXXnghfH3FihXMmDGDiy++mPT0dObPn096ejpPPfUUABkZGXzzm99k6tSpZGRkHPf86upqVq5cyV133cXo0aMZNWoUd911F4888gg1NTUAPP3005x33nkMGjQo/L3TTjuNOXPm8MlPfrLd2MeNG0d2dja//e1vA/rTiEzr/Emf0ZyczTx1JkXbi46rmCW6l3C073f07vxf5VNeXc6ysmXh+YA52TnsrtpNRW1ouMIwRmSM4GDNQXIeyQl/t3ltw4H9B1Ly5RJERFKueW5fZ+5N0jzA4uJiJkw4NtJSVlbGddddFz42M3JzcykrK4vpedu2baOmpoZJkyaFz02cOJHq6mq2b9/OhAkTKC0t5fTTT48r3pycHEpLS5kxY0Zc34+Fkj8JRE9Y4iTR5C4Z757wyISIi1Q3J3Jtm0ymnTSNrQe2RkwKn/iHJ4Ce8XchItIVnnjiCVasWMHvfve78LmqqiraTusaPnx4zEOtVVVVAK2e0fxz8zMOHjwYd4PM0KFDj5ujGDQN+0ogNAQZn/bmA667ah0QeSi6vfmJzYletOVmRES6RMGzrT+x3BuwoqIi5syZwzPPPMPEiRPD54cMGULbhs5Dhw4xdOjQmJ47ZMgQgFbPaP65+RkjRoyIe95eZWUlWVlZHd+YAFX+pNOat1aDjrdXg+O3d5NjHrvssYjb3TUncu1VDJuTwpbDyB1VEUVE+orCwkK+9a1vsWbNGs4///xW18466yxKS0vDx+7OK6+8EvMw6/jx48nIyKC0tJQLL7wQgM2bN5OZmcmpp54KQG5ubqt5hp3xl7/8pdWwdDIo+ZOEtDcv7cQTTkxxZD1HR/MBI4mUFF7z7DVRh4NFRFJm8KjIzR2DRwX+qiVLlnDnnXeybt26iM0Vc+bMYfr06Vx77bVMmTKFJUuWUFNTw5VXXhm+p7a2FnensbGR+vp6ampq6N+/P2lpaWRmZvLlL3+Z22+/nTPPPBOA22+/ndmzZ4cbRD73uc/x/e9/n+rqajIzM4FQkllbWxt+x9GjR6mpqSEtLY3+/fsD8MYbb1BeXs60adMC/3NpSYs8R6FFnmMT5PZqkhj9XYhIMvSkRZ7NjAEDBpCent7q/HvvvRf++dFHH+WOO+7gnXfeIScnh2XLlrVq4BgzZgy7du1q9f1rr72Wn/70p0Bonb+bbrqJJ598EoCrrrqKBx98MJzoAUyfPp0vfvGL4Srezp07Ofnkk4+Lt7CwMHzPd7/7XSorK3nwwQfb/f0FscizKn+SsHgqVx26b1z7/0qcvyPx5/dSSfm7EBHpQWIpas2ePZvZs2e3e33nzp1Rvz9o0CAefvhhHn744Xbvuf/++/nCF77ANddcQ1paGmPGjIka2/79+3nsscfoiqKTKn9RqPKXQndEWWD7Du28IiLSlQKr/EnCtL2b9D6Fl3W8PlRn1o8SERGRVpT8iYiIiPQhmvMn3UvzJOBow75JnCgsIiLS26nyJ9JHtN0/WERE+iYlf9I9tbf2UxLWhOortPOHiIiAhn2lu9JyLoHRzh8iItKSkj/pWlq/r8tF2oXl0pMv5dZP3priyESkr2neHlTbfqaWhn2la0VK/KKd78h940LNIW0/942LP8Ze5rHLHuOiky4CQokfwOCBg8nOzNY8QBGRJNmyZQvjx4+nrq4u5vs//vGPt9oCLlmU/EnyNa/dl4z1+4JOJtvTw5PMipoK8sfnU3RFEbPGz+JA9QFA8wBFpPe57bbbOOOMMxg6dCgf+tCHmDNnDhUVFa3uefTRRznllFMYNGgQ55xzDiUlx6bAVFdXM3PmTMaNG0e/fv24++67j3vHkSNHuP766xkxYgTDhw/nq1/9KtXV1a3uufXWW7nttttIS0sLv/O8885jxIgRZGdnc8kll/DnP/85fP8ZZ5xBbm5u1K3dgqLkT3qmrl4MuquSzCRZPHUxCyYvYHzWeBZMXsD/vvW/5DySQ3l1ORCaB5jzSA6TVk7q4EkiIvE72nCUrRVbkzra0L9/f1auXMmBAwcoKytj9+7dFBQUhK9v2LCBr33tayxbtoyDBw9y1VVXcemll1JZWQmE9gY+77zzWL58OWeffXbEd8ybN4+tW7eydetWtm/fzuuvv84tt9wSvr5t2zZeeuklvvjFL4bPVVVVceedd7J792727NnDxIkTufjii1sljddffz1Lly6lsbEx6D+WVpT8SfIVPHvs09F93VEv3FEkJzuHrPSs8LFhZGVkse6qdSmMSkR6u3fef4f36t5L6mjDPffcQ25uLmlpaYwcOZKbbrqJF154IXx9xYoVzJgxg4svvpj09HTmz59Peno6Tz31FAAZGRl885vfZOrUqWRkZBz3/OrqalauXMldd93F6NGjGTVqFHfddRePPPIINTU1ADz99NOcd955DBo0KPy9r3/963zmM59h8ODBpKen893vfpd3332XrVu3hu/51Kc+xbvvvssrr7ySpD+dEDV8yHHKj5Qzf/187v/0/WRnZqc6nMiCXgw6WiNK9qmdi60HeOyyx1i0cRFF24vC8wCnnTSt+/59i0iP0tzY0axkb0nEVQcMY9Lo0IhDsppAiouLmTBhQvi4rKyM6667LnxsZuTm5lJWVhbT87Zt20ZNTQ2TJh0bKZk4cSLV1dVs376dCRMmUFpayumnn95hXIMGDWLs2LHhc+np6YwbN47S0lImTpwY4++w85T8yXHyf5VPeXU5y8qWsXDywmAfPnhU+0lWMnXUZRxtWHf+juhJZizP74aa5wHOPHUmRduL1PQhIkkTadWBERkjOPGEE5P63ieeeIIVK1bwu9/9LnyuqqqKYcNa/3/68OHDw8O+HamqqgJo9Yzmn5ufcfDgQU477bR2n7F9+3b+8R//kR/96EcMGTKk1bWhQ4ceN0cxaEr++rB4/mUGCf7rLOhEKNZkMlpyF8Swbg+cE7h46uLwzwsmL0hhJCLS20T670Sk0YbACwwtFBUVMXfuXJ555plWVbQhQ4Zw+PDhVvceOnSIU045JabnNidrhw8fZvjw4eGfIZS4AYwYMaLdZPK1117jM5/5DLfeeis33njjcdcrKyvJysqK8M3gKPmTsFT9yywhHSWTQc3Xi5ZkJvKOHlgxFBGJR0VNBSMzRzIycyQ5I3OSOtpQWFjIt771LdasWcP555/f6tpZZ51FaWlp+NjdeeWVV5gxY0ZMzx4/fjwZGRmUlpZy4YUXArB582YyMzM59dTQNKHc3NxW8wyblZaWMn36dBYuXMjNN9983PWjR4+yY8cOcnNzY/2txkXJXx/W3r/MHt/+OGn906hrqEv6v8y6hYJnOx7WjZaIdZT8RUvw4qkYKmEUkR5o8dTF4RGnZI42LFmyhDvvvJN169bxyU9+8rjrc+bMYfr06Vx77bVMmTKFJUuWUFNTw5VXXhm+p7a2FnensbGR+vp6ampq6N+/P2lpaWRmZvLlL3+Z22+/nTPPPBOA22+/ndmzZ4cbRD73uc/x/e9/n+rqajIzMwF46aWXuPzyy/nBD37ADTfcEDH29evXM3r06KQnf+r2lVYqaiqYNX4WP7v0Z63Wg+uxYukyDuId0QQ9JNwDh5hFRLrKvHnzqKysZOrUqZxwwgnhT7MLLriAH//4x8yZM4dhw4axevVqnnvuufCQLYSqe5mZmbz44ovceeedZGZmMmfOnPD1Bx54gFNPPTX8GT9+PP/+7/8evn7aaadx7rnnsmrVqvC5BQsWcPjwYW655ZZWcb344ovhex5++GFuvvlm+vVLbnqmyl8EZnYFcEXLDpy+otPzwLqwCrWvsoabfr6ZB6/OZdSQ49vvE5KqRhQREQmUu3d4z+zZs5k9e3a713fu3Bn1+4MGDeLhhx/m4Ycfbvee+++/ny984Qtcc801pKWl8fzzz0d95muvvUZJSQmFhcnf+k7JXwTuvgZYk5eXN6fDm/u6LqxCLSnewZ92VrCk+A3u/vyZnftyR8ldoolqtOcnsnWdqnkiIj3SmWee2WoNv46cfvrpbNu2LYkRHaPkTzov1gaH5vviHHbNf2hj+Oej9Y1sfvsQAI+9vIstew4zcECoLL5q7rkdPyzZc+GiPb+j+YTtUeInIr1Mstbyk87RnD/pEfYcarFnorc57snaG1bul9a1cYiISJ+hyp90XstKXlA7bETQXNHbV1nDlB8emyvhQGV1HUuTMfcvGaINCbdXMSy8DHZt6Px7REREOqDkT1Kuo0aOJcU7aGwzgbfBPb65f6kQz5BzR8vP3HH4+HP3jYv8HS0BIyIiLWjYVxLTXrWpE1Wolo0ckZS+dYi6htbJX12DU7rrYMzv6BO0BIyIJFFjY2OqQ+jzgvo7UOVPEtPJilLLJg6IrZHjuXlTAgi0B9LyMyLSTQwePJg9e/YwevRo0tLSMLNUh9SnuDt1dXXs3buXwYMHJ/w8JX+SUpEaOU7OTvx/2L2ChmpFpJv4yEc+wv79+9m1axf19fWpDqdPGjBgAMOGDSM7OzvxZwUQj0jMWi7L0uMbOURE+oh+/foxatQoRo3SyENvoDl/kjLRGjl6gn2VNcx6aCP7qmpSHYqIiEjMlPxJyvT0Ro6OGlW6XADNNyIi0vtZLHvg9VV5eXm+adOmVIchKRJpCZrmhpWWjSpm8ImPDOepr58f9bs9UfmRcuavn8/9n76f7MzE55mIiEjymFmJu+d1dJ8qfyLtiFbZ62jHkW5VFWxe/6/t575xHX41/1f5lOwtYVnZsi4IVEREuoIqf1Go8td7RavqQeTK3sAB/Vg199xwo0pt/bH1ljIG9OO0Dw5l4IB+Ub+bEp1YLLpgbQEAJXtLcI7//wbDmDR6EqA9OkVEuhtV/kSi6KgyF62y116jSvM9HVUFe0KjSE52DlnpWeFjw8jKyGJC9oQURiUiIkHQUi8SuO463y3SfL2WC0vHspcwtN+okn1COku/lBv1u9A68exu29O1rOYt2riIx7c/Tlr/NOoa6ph20jQWTl6YwuhERCQISv4kcKlObjpKPjtaWLqjvYSj7Tiy4Kk/t/vdHXurOkw8u5OKmgpmjZ/FzFNnUrS9iP3V+1MdkoiIBEDJXy/W6U7N+8a1v51YO7tNdGa7NqBLkpzLl25gX1Xtcclny/l6zSJV5hJZgibad4dkDOhRO5osnro4/POCyQtSGImIiARJyV8v9pNXf0Lp3lKWlS2LbbguUuIX7XwEqUhu2jZq7KuqBSLvFdxRVQ8S20s42ndjSTyTQnsEi4hIC0r+erjm7syW2nZqrt62mtXbVrfq1AQofCf2pI7Cy1q89Nnwj91tu7aOks9ULiwdS+KZlPmS2iNYRERaUPLXC+Vk57C7ajcVtRVAqFNzRMYITjzhxKS+N5bkJhliadRoTqQSqeolKpbEs70ha+i+jTQiItKzKPnr4dpbay2uTs1o68G1qPa1J9XbtaUq+YxVtMQz/6GNEYesW+4akupGGhER6R2U/PVSqejUTGVVDVKffCYq0pB1tOVpWiaGnRZHc4+IiPQOSv56qbg6NXt4Y0Cqk89EtLc+4NH6RgYO6NfhwtGdFkBzj4iI9ExK/uQYVXxSpr0h6zM+PIxvXDg2YmK4r6om6tw/zREUEZFItL2bSDcQbcg62lxGaH+7uMuXbuCP/9f+FnYiItI3qfIn0g1EG7K+9IEXo85lbNsIEql55BsXjVX1T0REACV/It1ee4lh/kMbufI/Xoq4o0rbOYLqEBYRkWYa9pUOtTesKKkXqRFk6ZdyOVxd1/I0j296u/XfX3tNPD2kuUdEROKnyp90SOvLdU/tdQj/YO3Wjtc7VHOPiEifpeRPwlrukdss0vpyzXvlQuvt3aRrtdcI8r9b9/Xo9Q5FRCS5lPxJVB3tlSup016H8AeHZbL59otTFJWIiHR35m0qB3JMXl6eb9q0KdVhpEzzXrm19Y3hcxkD+rH+tqnqHBUREelmzKzE3fM6uk8NH9KujtaXExERkZ5HyZ+0q6fvlSsiIiLH05w/aVdP3itXREREIlPyJ9KHRN3v975x8P6+4780eJSWhhER6UU07CvSh0Td7zdS4hftvIiI9Eiq/In0AZH2+92y5zBPff38FEcmIiJdTZU/kT4i0pqNIiLS9yj5i8DMrjCz5YcPH051KCKBiLTfb2V1nfZrFhHpgzTsG4G7rwHW5OXlzUl1LHHRxH1pI9qajXd//kz2VdYwKkWxiYhI14o7+TOzHOBs4O+ADKAC2A783t21EFwqaeK+tNHRmo1Lincwz4cx0iJUuwcrLRQR6U06tb2bmX0M+BpwDTAaaAQOAbXAcGBQ07nfAf8JrHL3xshP6/56xPZuhZcdf27XhvbvP+mCYz8XPBt8PNKjNDeCbH77EABm8PJ3L9L2fSIiPVDg27uZ2X8CW4BPAIuAXCDD3Ue6+0fc/QRgFHAF8Gfgh8DrZnZBe88UkdRr2wii7ftERHq3mCt/ZvYgcJ+774rx/n7ATAB3XxV3hCnUIyp/kdwxLMo1NbHIMfsqa5jyw+eprT9WoM8Y0I/1t03tdPWv/Eg589fP5/5P3092ZnbQoYqISAcCr/y5+02xJn5N9ze6+6qemviJ9AXRGkE6K/9X+ZTsLWFZ2bKgwhMRkSRQt29vNHhU+92+Ii101AgSTcHaAgBK9pbgHHvG6m2rWb1tNYYxafQkCqcXBhu0iIgkJKHkz8yWu/sNQQUjAWlnOZd9lTXc9NDGyPu6Sp/03LwpCT8jJzuH3VW7qaitAMAwRmSM4MQTTkz42SIiErxEK3/TA4lCusSS4h38aWdFeG03kUS0rOgt2riIou1FGAbAtJOmsXDywlSFJiIiUXSY/JlZQ3uXgNjXiZEuk//QxuPOtVzOo3lf14EDWk/5XDX33C6JT3qfipoK8sfnM/PUmRRtL2J/9f5UhyQiIu2IpfL3DpDr7uVtL5jZ28GHJMkQaV/Xk7MHpy4g6VUWT10c/nnB5AUpjERERDoSS/L3DDAeOC75A9YGG44EoW0Fr3k5j2bN+7ou1dw/ERGRPqfDpV7c/Z/cPeKWEe7eM/e+7WOCXM5DREREeraY1/mTniuR5TxERESkd+l0t6+ZrVDFr2cCg07AAAAgAElEQVQJYjkP6cPuG9f+upHtLCskIiLdVzyVv4sCj0JEuq9IiV+08yIi0q1p2FdERESkD1HyJyIR7ausYVaENSNFRKRnU/InIhE17wgjIiK9i5K/bqz8SDnXrb1OuyVIl8p/aCNX/sdLrHz5LVx7+IiI9DrxJH97Ao9CIsr/VT4le0tYVrYs1aFIH9NyR5hyHxb5psGjuigaEREJkrn+ad+uvLw837RpU1LfUbC24LhzJXtL8AjbJhvGpNGTWp0rnF6YtNikb2reEaa2vjF8LmNAP9bfNlU7woiIdGNmVuLueR3dF3Plz8w0RNxFcrJzyErPCh8bRlZGFhOyJ6QwKukrtCOMiEjv1plFnl83swXuXpS0aPqg9ip3izYu4vHtj5PWP426hjqmnTSNhZMXdnF00hdpRxgRkd6tM8nfHmCVmf0L8F13X5ekmASoqKlg1vhZzDx1JkXbi9T0IV1GO8KIiPRunZrzZ2Z/D9wNnAe8CHzH3X+fnNBSryvm/ImIiIgEIfA5fwDu/oK7XwBcCmQCL5rZM2aWE2ecIiIiItKF4mricPe17n42MAM4EdhsZivN7GOBRiciIiIigUqog9fdf+nuucDVwERCTSE/NrMPBhKdiIiIiAQqkOVb3H01cCbwPeA6YEcQzxURERGRYHWm2xcAMxsAjAXGAx9v8+twwICjAcYoIiIiIgGJOfkzs18SSvJOBvoTSvKqgK3A68Avm359HXgz8EhFREREJGGdqfwNBYo5luC97u7a51ekj9lXWcNNP9/Mg1fnars3EZEeKObkz92nJjMQEekZLl+6gX1VtSwpfoO7P39mqsMREZFO6vScv2ZmNhL4AvBRoBJY7u4HggpMRLqX/Ic2crS+kX1VtQA89vIutuw5zFNfPz/FkYmISGfElfyZ2fnAc8D7wBvA+cA64ICZfQXY6O7aBV6kl9lzqPrYgbc5FhGRHiHepV4WA78mVPW7kFDzR7MzgNsTjEtEupmlX8rlcHVd+NiByuo69lXVpC4oERHptHiTvzOAFe5eT+i/AS39kdDevyLSiywp3kFjm73AG9xZUqwiv4hITxJv8reH0JIvkZQD2uFDpJcpfesQdQ2tk7+6Bqd018G4n1l+pJzr1l7H/ur9iYYnIiIxirfhYzlwh5k9D+xsc20MEP9/DUSkW3pu3pTIF+4bB3fsO/784FEwP/pmP/m/yqe8upxlZctYOHlhAFGKiEhH4k3+/g34LFAKPEpo6HeomeUS2uLt+WDCE5Fu7/0IiV+L8wVrC467VLK3BG8xY2T1ttWs3rYaw5g0elL4fOH0wmBjFRGRuId9Gwklf/cAswg1fBQDm4AjwG2BRCfH2VdZw6yHNmqSvaRe4WWhT0f3RJCTnUNWelb42DCyMrKYkD0hyAhFRCSCeCt/P3L3W4AfmNkPgRxC8/zeAba4e0NQAUprS4p38KedFVpgV3qM9qp3izYuomh7Eda0WMC0k6Zp6FdEpAvEm/zNNrMt7v5f7u7Aq00fAMxstrs/GkiEfVD+Qxsjnj9a38jmtw8BxxbYHTjg+OLtqrnnJjU+EQAKng39esewju+JoKKmgvzx+cw8dSZF24vU9CEi0kXiTf5mAc+Y2Q53X9980sz6EVoD8EZCcwElDg/t/RLDG4/vmTnAcCbx49BB0wK7J2cP7uLoRIKxeOri8M8LJi9IYSQiIn1LXMmfu/+vmf0L8ISZne3u/2dmHwAeB04nNB9Q4hQp8QP4AIfCPzcvsLv06lxGDcnooshEIhg8KnLTx+BRXR+LiIh0KObkz8z6t5zL5+4PmtkE4Fkz+zrwU0Jr/H3S3d8KPNLeqqMJ81E01B/V3D9JvQ6WcxERke6lM92+75tZiZn9p5nd1LS/778QSvh+S2h5lwuU+HWdOgYktMCuSDKoI11EpHvrzLDvdcAE4BPAJYS6exsJLeh8EHgT+KyZ/dnd/xpwnL1XpAnxUSbQ77w3/kqhSFe4fOkG9lXVqiotItJNxZz8ufsvgF80H5tZNqFE8Kymz0xgITDAzN5z96EBxyoi3VDL7vSj9Y3sq6oFju9IVxe6iEj3EG+3L+6+n9Bw72+bz5lZGnAmoQqhxEsT6KWH2nOo+tiBOtJFRLolCy3TJ5Hk5eX5pk2bUh2GSI+wr7KGKT98ntr6xvC5jAH9WH/bVHWki4h0ATMrcfe8ju6Ld3s3EZFWlhTvoLHNPyYb3FlS/EaKIhIRkUiU/IlIIErfOkRdQ+vkr67B1ZEuItLNxD3nT0SkpefmTUl1CCIiEgNV/kRERET6ECV/IiIiIn1IQsmfmS0PKhARERERSb5EK3/TA4lCRERERLpEhw0fZtbQ3iVAiwSKiIiI9CCxVP7eAf7O3fu3+fQD/pbk+EREREQkQLEkf88A49u5tjbAWEREREQkyToc9nX3f4pybU6w4SSHmWUDa4CjwCDgu+7+P6mNSkRERKTr9ZVFng8CF7h7g5mNBX4OKPkTERGRPqfT3b5mtiIZgSSTuze4e3PjyhCgLJXxiEiwyo+Uc93a69hfvT/VoYiIdHvxLPVyUeBRAGb2RTN70cwqzaw+wvX+ZnafmZWbWZWZPdE0nBvr8082sw3Ab4Cng4xdRDpvX2UNsx7ayL6qmoSflf+rfEr2lrCsbFkAkYmI9G7m3rnVWszsr+7+scADMfsskAVkAsvdfUCb698DriW0tuAB4GFgkLtfYmYnAC9EeOwad7+zzXPGAC+4+5iOYsrLy/NNmzZ1+vciIh1b8NSfeeyPb3HNOSdx9+fPjHhPwdqCqM8o2VuCR1hxyjAmjZ7U7vcKpxd2LlgRkR7AzErcPa+j+7rNnD93XwdgZn/fzi03AIvc/a9N930beMPMxrj7TqDd36yZpbt7bdNhJVAVVNwiEpv8hzaGfz5a38jmtw8B8NjLu9iy5zADB4QGIlbNPTfmZ+Zk57C7ajcVtRVAKOkbkTGCE084McDIRUR6l26T/EVjZsOAjwIlzefc/U0zqwQmADs7eMREM/sB0EDo9/zPUd51A6FEk49+9KOJBS4iEe05VH3swEPHJ2cPPu6+WCp0izYu4vHtj5PWP426hjqmnTSNhZMXBhmuiEiv0iOSP2Bo06+H25w/1OJau9x9I/CpWF7k7suB5RAa9u1EjCISRXNFb19lDVN++Hz4vAOV1XUsvTqXUUMyOv3cipoKZo2fxcxTZ1K0vUhNHyIiHYgn+dsTeBQdax6mHdbm/HBCw7gi0kMsKd5BY5u5xg3uLCl+o925f9Esnro4/POCyQsSjk9EpLfrdLevu09JRiAdvPMQ8BYwsfmcmX2MUNXv1a6OR0TiV/rWIeoaWid/dQ1O6a6DKYpIRKRv6TbDvmbWH0gDBjYdN4//1HqoJXk5cJuZPU+o2/cHwLqmZg8R6SGem9fl/34UEZEW4lnnDzO7yMw+HXAsXwGqgXVA/6afq4GTmq7fS2iLtj8RGnruD3w54BhEREREerW4kj9Cidhnmw/M7HQz+72ZvWxmX4nnge7+U3e3CJ+dTdcb3P1Wd8929yHuPsPdNbNbREREpBPiTf5OBTa0OP43YDShpVhWmNnliQYmIiIiIsGLd85fI3AEwMxGA9OA6e7+WzN7F/gO8KtgQhQRERGRoMSb/P0FuIDQlmpfJLTeXnHTtfXALQlHJiK9133j4P19x58fPArm7+j6eERE+pB4h33/FfiemRUBdwE/82ObBH+AUGVQRCSySIlfO+f3VdYw66GN7KuqSXJQIiJ9Q1yVP3dfa2aXAFcDhUDLlVWnAtsDiC1lzOwK4IqxY8emOhSRnq/wsvjvL3iWy5duYF9VbdyLQIuISGtxr/Pn7i8QGvaN9MxV8T63O3D3NcCavLy8OamORaSv2fLOsV0cF/zHS+yrqgXgsZd3sWXPYQYO6BfeKk5ERDov8EWe3f3GoJ8pIj1YwbPHn7uj7U6Nxyz6wH3hn/fsf//YBYc9h6o5OXtwkNGJiPQ58c75ExFJilVzz2XV3HNZ+qVcDlfXhc87UFldx9Krc9v9ruYHioh0TMmfiHS9waM6PL+keAeN3noP4AZ3lhS/0e5jlxTv4E87K6LeIyLS13WbvX1FpA+JYTmX0rcOUdfQOvmra3BKdx0k/6GNx91/tL6RzW8fAlrPD2ymeYIiIiFK/kSkW3pu3pR2r0VK/vYcqj52kOT5geVHypm/fj73f/p+sjOzk/IOEZFkUfInIj1O2yrevsoapvzw+fBxy/mBo4ZkBP7+/F/lU15dzrKyZSycvDDw54uIJFNMyZ+ZXQmcCKxz920tzt/k7g8mKzgRkVhEmx8YaW3AgrUFcb2nZG8JzrH3rN62mtXbVmMYk0ZP6tSzCqcXxhWDiEiiOmz4MLN7gXnAWOB/zOyfW1y+PlmBiYjEKtr8wGhqazMo3fxpamvTY3pPTnYOWelZ4WPDyMrIYkL2hM4HLSKSIuZt/rV83A1mfwZy3b3ezD4AFAEl7j7fzDa7e/vrLvRweXl5vmnTplSHISJJcva//pZ9VbV8efJJMe8esmjjIoq2F2EYADPHz9TQr4h0C2ZW4u55Hd0Xy7BvP3evB3D3A2Y2HXjMzP4LLRUjIj1E2yaRo/WNEXcPadZed3BFTQX54/OZeepMirYXsb96f/KCFhFJgliSv3fMbKK7lwK4+1EzyweWA71yo03t7SvS+8XbHbx46uLwzwsmL4hyp4hI9xTLsO9HgHp3fzfCtfPc/ffJCi7VNOwr0js1dwfX1jeGz2UM6Mf626YmpTtYRKQrxDrs2+GwrbvvjpT4NV3rtYmfiPRe8eweIiLSW3R6zp6ZrUhGICIiXSXe7mARkd4gnkWeLwo8ChGRLhRt9xARkd5O3boiIiIifYiSPxEREZE+RMmfiIiISB+i5E9ERESkD4mn4WNP4FGIiHTGfePg/X3Hnx88Cubv6Pp4RER6kE5X/txdbXIiklqREr9o50VEJCyeyp+ISNcovCyY7xQ8m3gsbeyrrOGmn2/mwatztSuIiPQocSV/ZpYOFADjgQrgL8Cr7v5mgLGJiHRbly/dwL6qWpYUv8Hdn++V25yLSC/V4d6+Eb9k9gTweUJJ32BgDGDA+8AWoMzdbwwuzK5lZlcAV4wdO3bOjh2aPyTS7dwxLMq1w4G8Iv+hje1eO1rfyOa3DwFgBp/4yHAGDog8i2bV3HMDiUdEpCOx7u0b77DvxcDN7v7jppdlAjnAhBafHsvd1wBr8vLy5qQ6FhFJsnaaRx7qN4K5o38e8St7DlUfO/DQ8cnZg5MVoYhIoOJN/t4C/q/5wN2rgT82fUREkmvwqPa7fTurnSaR4Y0HI1bt9lXWMOWHz4ePHaisrmOp5v6JSA8Rb/J3L/BPwK8DjEVEJDadXc4lnsaRdr63ZPi9NLaZLtPgntDcPzWPiEhXimuRZ3f/b2Cnmf2PmV1oZmkBxyUi0i2VvnWIuobWyV9dg1O662Dcz7x86Qb++H8VLCl+I9HwREQ6FG/Dx7eA+1qcqgO2AmVNn1fd/X8CiTCF8vLyfNOmTakOQ0SSKcnNI9EaR6BzzSOQnAaS8iPlzF8/n/s/fT/ZmdmBP19EukasDR/xbu/2PWAloS7fM4DZwLNAFjAPWBvnc0VE+pRIzSNdLf9X+ZTsLWFZ2bIuf7eIdL14K397gS+5+/+2c324ux9KNLhUU+VPpA9I0VZxBWsLqK3N4A8vX0JjY//w+X796pl8zq9JT69N2rublewtwTn+vwGGMWn0pKS/v3B6YdLfIdKXJHupl5XARUDE5K83JH4i0kekcC/gnTtPo+2/v92NnbtOY/ypryT9/TnZOeyu2k1FbQUQSvpGZIzgxBNOTPq7RSR14k3+dgH/bGZ/A37i7g0BxiQi0usVTi/k0m0v8jevbHXevT/DfSKF0+d1SRyLNi6iaHsRhgEw7aRpLJy8sEveLSKpEW/ydw8wCFgKLDKzDcArNDV8aJs3EZGOPTdvSqpDoKKmgvzx+cw8dSZF24vYX70/1SGJSJLFO+fPgI8R2skjh2O7e5xCqInkfXcfEmCcKaE5fyIiItJTJHXOn4cyxjebPk+1eGkGcGbTR0RERES6mbiSPzP7byDN3b/Y8ry71wCbmj4iIiIi0s3Eu87fRbSztZuZfd/MZsUfkoiIiIgkS7zJ3wjg7Xau7Qb+Jc7ndgtmdoWZLT98OPHV/UVERES6k3iTv+3AxHauvQaMi/O53YK7r3H3G4YNi7Ltk4hID7WvsoZZD21kX1VNqkMRkRSIN/n7KfAdMzs1wrUPAUfijkhERJLq8qUb+OP/VbCk+I1UhyIiKRDvOn8PAJ8CNpnZUuBp4B3gNOBOYH0w4YmISEv5D21M6PtH6xvZVxXaOu6xl3exZc9hBg7ofB1g1dxzE4pDRFIn3qVeGs1sBnALMJ9jc/wM2ALcGkx4IiI9XIr2Dm7PnkPVxw48dHxy9uAuj0NEUifeyl/zWn8/MrN/I7TA8weBd4E/a7s3EZEmkRK/aOc7kEjFbV9lDVN++Hz42IHK6jqWXp3LqCEZcT9XRHqWuJO/Zk1JYFnTR0Skdyu8rGufVfBsYK9bUryDxja7OjW4s6T4De7+vNbmF+kr4m34EBGRHqb0rUPUNbRO/uoanNJdB1MUkYikQlx7+/YV2ttXRBJ2R5Qlo+7QWqIiEpxY9/ZV5U9ERESkD1HyJyKSTINHde68iEiSxdXwYWYXAfXu/ruA4xER6V1SsJyLiEg08Vb+7gU+23xgZqeb2e/N7GUz+0owoYmIiEhXy39oY8KLiUv3Fm/ydyqwocXxvwGjgRJghZldnmhgIiIikJy9iK959hrO/dm57K/eH9gzU61gbQEFawsSfs7R+kZe+1ul9n7uxeJN/hpp2r/XzEYD04C57v5PwD3Ad4IJT0RE+rpk7EX8zvvv8F7deywrWxbYM3uLPYeqqaqt197PvVhcS72Y2YvAOne/28zmAQuBke7uZvb3wNPuPjzYULuelnoREYksiGHBnQPv7/CehvpM9v91JvgAsHqyP1ZE/wHVHX6vPUdsB9jx/90zjEmjJ8X93I4UTi9MynNb/j3838B7qbV3GFu7iAFEWWIoiqP1jWx++xAAZvCJjwyPa+/ntrQXdNdI9lIv/wp8z8yKgLuAn/mxLPIDhCqDIiIiCXl//ydC+9ABeNNxAjL9ZPr7kBbPNPr7ECZkT0joud3BUTtAI9WUD/hV3M+ItPez9D5xL/LcVOG7GqgGFrp7ZdP5B4E8d58cVJCposqfiEjqNO9FXFt/rJ6QMaAf62+bmtBexIs2LqJoexGGATBz/EwWTl6YcLzRBDEXrz0le0twEq9m1tZm8IeXL6GxsX/4XL9+9Uw+59ekp9cmrXopwUn6Is/u/oK73+Du85oTvyYDgFXxPldERASi70WciIqaCkZmjuS0rNOYNX4WB6oPJPS8VMvJziErPSt8bBhZGVmdrmbu3HkabetB7sbOXacFEaZ0I3Gt8wdgZp8DxgH7gS3AX9y92t1vDCq4VDGzK4Arxo4dm+pQRET6rGTtRbx46uJwJW7B5AUJPStWya6ata1mTjtpWqermZdue5G/tarlgHt/hvtECqfPCyxWSb14Gz6WA18F3gWGA5lAA/BX4FWgzN3vDjDOlNCwr4iIBGFfZQ03/XwzD16dm9CQdXv++fl/5tXyVxmZOZKckTnsr97P4qmLA3+PdG/JHvadBdzu7h9298GEKoD5wM+B/sC1cT5XRESk10nGcjUtLZ66mJOGnsSgtEEsmLxAiZ9EFe+wbxXwh+YDd38TeBN4MoigREREulIyd7Q4Wt/IvqpaAB57eRdb9hwOZPkU0BIqEp94k7+fAtOB4uBCERER6X0iLZ9ycvbgwN+jblyJVbzJ327gW2b2NvBjd68PMCYREZEulawKWvNyNc0cqKyuY2mS5v6JxCLeuvOPgLHAYqDczJ4xs7vM7CozU4usiIgIyVuuRiQR8Vb+hgAfAyYAOU2fWYT29O1nZu+7+5BgQhQREemZkrVcjUgi4kr+mrZya27yeKr5vJllAGc2fURERPq05+ZNSXUIIseJe5HnSNy9BtjU9BERkZ7qvnHw/r7jzw8eBfN3dH08IhKYRHb4GAZcAnyY0GLPG9x9V1CBiYhICkVK/KKdF5EeI67kz8wmAL8BRgKVwDDAzezXwFx33xNciCIifVjhZamO4HjJjKng2eQ9W0SA+Lt9lwKbgZHuPgI4AfgHQsngH8zsgwHFJyIiIiIBindv3yrgSnf/bZvz/YF1wB537/FbvGlvXxHps+4YFuXa4a6LQ0Riluy9fSuAD7Q96e4NhNb+uyTO54qIiIhIEsXb8LEKuN3MfuPubRcrsgSeKyIi3cHgUe13+0r3oa5siUO8Sdr/Az4F/MXMHgR+S6jj92PAXcCLwYQnIiIpocShZ1BXtsQh3uTvHODvgUXAt4G7m84bUArcnHBkIiIiPV0qu7WDerc6sHudeJO/YuBcd/+2mX0POAvIBt529y2BRSciIiIigYo5+TOzAe5e33zYfN7d62ixo4eZnQM85e4fCixKERGRnijZVbNoXdmq2Ek7OtPt+z0ze9fMfgs4cJmZnW1mg9rcN5AIncAiIiIiknqdGfb9GXAQmECo8jcfWAg0mtlfgTJgK6FGkO0BxykiIiJtqSu7SxWsLQCgcHphiiNJTMzJn7vvAHYAmNk/AFcQ6vCd0OLzD4S2e7sx8Ei7kJldAVwxduzYVIciIiLSPnVlSxziavhw95b/pNgFrAkmnO7B3dcAa/Ly8uakOhYREQmY1saTPi7eHT5ERER6Jq2NJ31c3DtxmNkwQtu4fZjQ8O8Gd98VVGAiItIHpXJdvGS8Xx230g3FlfyZ2QTgN8BIQnP8hgFuZr8G5rr7nuBCFBEREZGgxFv5WwpsBq5x94qm5V6mArcDfzCzs939naCCFBGRPqIrKmVaG0/6uHjn/E0EfuTuFQDufsTdnwXOA7YB9wYUn4iIiIgEKN7kr4IICzm7ewOwmNBcQBERke6nvTXwtDae9BHxDvuuAm43s9+4+8E21yyB54qIiCSXlnORPi7eyt//A6qAv5jZd8zsk2Z2opl9GrgLeDGwCEVEREQkMHElf+5eDfw98BjwbeAPwE7geaAeuCmY8EREREQkSHEPz7p7DfBtM/secBaQDbzt7luCCk5EREREgpXIIs+fA8YB5cBrwO+aKoIiIiIi0k3Fu8jzcuCrhHb2GA5kAg1m9lfgVaDM3e8OLEoRERERCUS8DR+zgNvd/cPuPphQBTAf+DnQH7g2oPhEREREJEDxDvu+R6jJAwB3fxN4E3gyiKBEREREJDlirvyZWf8Wh48AlwYfjoiIiIgkU2eGfd83sz81zffbD1xpZvPaJIUiIiIi0o11Ztj3OmAC8AngMuCDwL8T2unjdzQ1egCvNg0Di4iIiEg3E3Py5+6/AH7RfGxm2YQSwbOaPjOA7wIDzOx9dx8ScKwiIiIikqBEFnneD/y26QOAmaUBZxKqEIqIiIhIN9OZho+RZrbKzCrMrNrMNprZFS3vcfc6d9/s7o8EH6qIiIiIJKozDR8PAJ8H1gKPAoOBp83sS8kITERERESC15lh32nAt9z9weYTZrYM+D6hxZ1FREREpJvrTOUvG/hTm3P3Aiea2UnBhSQiIiIiydLZ7d0a2xzvBgwYEUw4IiIiIpJMne32/Xcz+yOwpemzvem8BRqViIiIiCRFZ5K/Bwgt43I1MArwFtfuM7PfE1ro+VVgh7v78Y8QERERkVTqzCLP32z+2cxGElrLbwKQ0/T5FpBJKCmsBk4INFIRERERSVhcizy7ezlQ3PQBwMwMGEcoITwzkOhSpGn9wivGjh2b6lBERESkG8p/aCMAq+aem+JIOq+zDR/t8pDt7v64u98R1HNTwd3XuPsNw4YNS3UoIiIiIoHqzA4fXzGz/p15uJmNNbMpnQ9LRERERJKhM5W/bwFvmtldZnZWezeZ2QfM7BozWwNsBj6YaJAiIiIiEozONHx8wszygZuB75nZe8DrwH6gFhgOnAx8FDgIrARudPc9gUctIiIiInHpVMOHu68CVpnZKYS2e5sI/B2hfX73AuuBl4AX3L0u4FhFREREJEHxdvu+CbwZcCwiIiIikmSBdfuKiIiISPcXU/JnZlea2TfMbHyb8zclJywRERERSYYOkz8zuxeYB4wF/sfM/rnF5euTFZiIiIiIBC+Wyt9lwDR3/waQC/yDmd3XdM2SFpmIiIhId/Pun6HwslRHkZBYkr9+7l4P4O4HgOnAGDP7rxi/LyIiIiLdRCzJ2ztmNrH5wN2PAvmA08P38BURERHpa2JJ/q4D/tbyhLs3uvs/Atq6TURERKQH6XCdP3ffDWBm/dy9sc213ycrMBEREZHu5GjDUf7KUfbTkOpQEtKZOXuvm9nMpEUiIiIi0o298/47vIezzA6nOpSEdGaHjz2Etnb7F+C77r4uSTGJiIiIdJmCtQVRr5fsLcHx0IHBat6DjBvAjYK1k9r9XuH0wiDDDEzMlT93vxC4EKgGfm1mvzOz85IWmYiIiEg3kJOdQ1Z6VvjYHPr7EDL95BRGFb9O7e3r7i8AF5jZdGAR8KKZPQt8z93/nIT4RERERJIqlgrdoo2LKNpehDUVAIc2TOSD9ddQOP3cJEcXvLjW6XP3te5+NjADOBHYbGYrzexjgUYnIiIi0g1U1FQwMnMkp5HGLE6g3ipTHVLcElqk2d1/6e65wNXAREJNIT82sw8GEp2IiIhIN7B46mJOGnoSg+jHAs/ixLqvpTqkuAWyQ4e7rya04PP3CK0LuCOI54qIiIhIsDo15w/AzAYAY4HxwMfb/Dqc0H6/RwOMUUREREQCEnPyZ2a/JJTknQz0J5TkVQFbgdeBXzb9+jrwZuCRioiIiEjCOlP5GwoUcyzBe93d9yQlKhERERFJipiTP3OkIf4AABJRSURBVHefmsxARERERCT5Oj3nr5mZjQS+AHwUqASWu/uBoAITERERkeDFlfyZ2fnAc8D7wBvA+cA64ICZfQXY6O5vBBaliIiIiAQi3qVeFgO/JlT1u5BQ80ezM4DbE4xLRERERJIg3uTvDGCFu9dD807HYX8EtOeviIiI9FpH6xt57W+V7KuqSXUonRZv8reH0JIvkZQD2uFDREREeq09h6qpqq1nSXHPm+UWb8PHcuAOM3se2Nnm2hjgYAIxiYiIiKRc/kMbWx3vHFjJSXX1vLivnn3v1QLw2Mu72LLnMAMHHF9PWzX33C6Js7Pirfz9G6HFnUsJzf9zYKiZ5RLa4u35YMITERER6V5+Xnv+sQMPVQF7kngrf43AZ4FbgVsINXwUN10rA25LPDQRERGR1GlbuStYO5Ta3cMprj8rfM6Byuo6ll6dy6ghGV0cYXzirfz9yN0b3P0HwN8BnwAuBXKBT7r734IKUERERPqmgrUFFKwtSHUYrezc/2ka27S6Nrj3qLl/8SZ/s83sqwAe8qq7r2v6tcHMZgcYo4iIiEi3UFl9InVtBk7rGpzSXT2n3SHeYd9ZwDNmtsPd1zefNLN+hOYA3gg8GkB8IiIiIt3GJ0/+CYU+mvyjC4Du29QRTVyVP3f/X+BfgCfM7GQAM/sAoXl/+YTmA4qIiIhINxNz8mdm/Vseu/uDwFPAs2Y29f+3d+/BdpXlHce/jwkQBBJSEkTlfvEPBDpiKDBqC5VKEJGWqci0WIhTUzqKOiiXP0ADOlMs1lLo0CFTOBRaCWWQShpLIIxQEERuym2AEAjXCRASkhBCQsLTP9ba4bC7zzn75OzrWd/PzJ6d9a53v++zJ5M1v7zrsinu/N2B4po/7/aVJEnqQaM57bs2Ih4DHgJ+W76fQxEAFwHXAKdlZv896lqSJKkiRhP+TgUOoriz9xiKX/F4l+KBziuBJcDREfFIZj7T4jo7KiKOA47bd999u12KJEkaQe1hzG29/u6i/RhY+2q5sZjrmFm27wxnLm7fvG3QdPjLzHnAvNp2REyjCIK/X76+BJwHTIyINzNzcotr7ZjMnA/MnzFjxte6XYskSeoBm4Nfk+09bEvv9iUzl1Oc7l1Ua4uIrYADKFYIJUmSet/Ase39/KwFYxu/xbY4/DWSme9QXAv4UCvHlSRJUmu0NPxJkiT1nWZW5uZMGdvne8iW/sKHJEmS+pDhT5IkaSTb7Ty69h7maV9JkqSRnLmYWTfPgmWPVPPn3SRJktSfDH+SJEkVYviTJEmqEMOfJElShRj+JElSX9uw8V0ef3k1r655u9ul9AXDnyRJ6kkbNm3giRVPsHzd8mH7vfTGOtas38gltz3docr6m496kSRJbfXly+/Zos89uc1zbOJNTvjpHD688S8b9tmw8V1eXbMegP+49zkee2kVW09szdpWPz7GpRmGP0mSNKJZN8/a4s8u3Xr1qPq/FYshcvP2yol3sHLiHZDBB3O/9/VdveJwYD9gIpkbeWL1PUzeZXRhc88N3x1V/35n+JMkSW21/4cnj6r/hk0H8uKaF1mxfgUAQTB10lR22343tpqw1eZ+69dP4teLPwZMKFsmsn7Nfuyz/9Nss836pucbmDk+V/iGYviTJEkjGpg50NH5LrjnAq5/6nqCAOCoPY7ivMPOe1+fc298hPviBd7lvVXCCbEVe+c3+OHMAzpabz8x/EmSpJ6z4u0VTN92OtO3nc6B0w9seNPHg8+/wTub8n1t72xKHnxuZafK7EuGP0mS1HMuPvLizdcZnnvYuQ37/OJbnwHeu6FkvN6g0Wo+6kWSJKlCDH+SJEkVYviTJEmqEMOfJElShRj+JEmSKsTwJ0mSVCGGP0mSpAox/EmSJFWI4U+SJKlCDH+SJEkVYviTJEmqEMOfJElShRj+JEmSKmRitwuQJEkai+v+5vBul9BXXPmTJEmqEMOfJElShRj+JEmSKsTwJ0mSVCGGP0mSpAox/EmSJFWI4U+SJKlCDH+SJEkV4kOeJUmSRqmfHyztyp8kSeo9A8fCske6XcW45MqfJEnqSQP5IZg50O0yxh1X/iRJkirE8CdJklQhhj9JkqQKMfxJkiRViOFPkiSpQgx/kiRJFeKjXiRJkpq1y4F9//gZV/4kSZIqxPAnSZJUIYY/SZKkCjH8SZIkVYjhT5IkqUIMf5IkSRXio14kSVLvuGg/WPvqe9tzphTv2+0MZy7uTk3jjCt/kiSpdwwOfs20a9Rc+ZMkSUMbOLbbFbynXbXMWtCecXuUK3+SJEkV4sqfJEkaWqdXxWrX+DVSsRW6dnHlT5IkqUIqFf4iYqeIWBkRJ3e7FkmS1MB2O4+uXaNWtdO+5wJ3dbsISZI0hNrjXGo3d3iqt+Uqs/IXEfsAOwEPdLsWSZKkbumZ8BcRJ0XEnRGxOiI2Ntg/ISIuiojXImJNRNwQEdNGMcUPgPNbV7EkSVL/6ZnwB6wELgO+PcT+c4DjgUOBXcu2awAiYvuIuL/B6/vl/sOB1zNzSXu/giRJUm/rmWv+MnMhQEQcMUSX2cAFmflM2e8s4OmI2DMzlwIzhhn+EOCgiLgZ2BdYGxFLMvOeVtUvSZLUD3om/A0nIqYAuzPoer3MXBIRq4GDgKXDfT4zLwEuKceaAzw9VPCLiNkUQZPdd9+9BdVLkiT1jr4If8Dk8n1VXfsbg/Y1JTPnjLB/LjAXYMaMGTmasSVJknpdL13zN5w15Xv9Y793BFZ3uBZJkqS+1Rcrf5n5RkQ8DxwM/BYgIvamWPV7uJu1SZKkahiYOdDtElqiZ1b+yke5TAK2Lrcnla8ou8wFzo6IvSJiMvAjYGF5s4ckSZKa0DPhD/gKsA5YCEwo/7wO2KPcfyEwH7gPeKns48+0SZIkjULPhL/MvCozo8Frabl/U2Z+NzOnZeYOmXlCZi7vctmSJEl9pWfCnyRJktrP8CdJklQhhj9JkqQKMfw1EBHHRcTcVavqnyktSZLU3wx/DWTm/MycPWVK/TOlJUmS+pvhT5IkqUL64hc+JElSxcxa0O0Kxi1X/iRJkirE8CdJklQhhj9JkqQKMfxJkiRViOFPkiSpQgx/kiRJFWL4kyRJqhDDnyRJUoUY/hrwt30lSdJ4ZfhrwN/2lSRJ45XhT5IkqUIMf5IkSRVi+JMkSaoQw58kSVKFGP4kSZIqxPAnSZJUIYY/SZKkCjH8SZIkVYjhT5IkqUIiM7tdQ8+KiNeA59o0/BSgH34/rtt1dnL+ds3V6nFbMd5YxpgGLB/j/Gqtbv877aR++K69UGOnamjnPK0cu9vHTejMsXOPzJw+UifDX5dExNzMnN3tOkbS7To7OX+75mr1uK0YbyxjRMT9mTljLPOrtbr977ST+uG79kKNnaqhnfO0cuxuHzfLz/fMsdPTvt0zv9sFNKnbdXZy/nbN1epxWzFet/9e1VpV+vvsh+/aCzV2qoZ2ztPKsT1uDuLKn6RR6aX/vUpSv+ilY6crf5JGa263C5CkPtQzx05X/iRJkirElT9JkqQKMfxJkiRViOFPUktExMkRcU/5+my365GkfhARiyJieUSc26k5J3ZqIknjV0TsCJwBHA5sD/wyIj6RmZu6W5kk9bxTgaOAXTs1oSt/klrhUODOzFyfma8DS4F9uluSJPW+zHyx03Ma/iQBEBEnRcSdEbE6IjY22D8hIi6KiNciYk1E3BAR08rdOwErB3VfWbZJ0rg2xmNnVxj+JNWsBC4Dvj3E/nOA4ylW+WqnJ64p318Hpg7qO7Vsk6TxbizHzq7wmj9JAGTmQoCIOGKILrOBCzLzmbLfWcDTEbEncC/wdxGxDbAdsBewpM0lS1LXjeXYmZlLO1JkHcOfpBFFxBRgd+CBWltmLomI1cBBmXlTRFwM3F7uPsObPSRV3UjHTmBpRFxJsSq4TUQcmpnHtbsuw5+kZkwu31fVtb9R25eZVwNXd7IoSepxzRw7v9rRivCaP0nNWVO+T6lr3xFY3eFaJKlf9OSx0/AnaUSZ+QbwPHBwrS0i9qb4n+vD3apLknpZrx47DX+SgM2PI5gEbF1uTypfUXaZC5wdEXtFxGTgR8DCbl2wLEm9oB+PnYY/STVfAdYBC4EJ5Z/XAXuU+y8E5gP3AS+VfU7ufJmS1FP67tgZmdnN+SVJktRBrvxJkiRViOFPkiSpQgx/kiRJFWL4kyRJqhDDnyRJUoUY/iRJkirE8CdJklQhhj9JfSci5kREDvEaVw+ejogPRsSyiPijQW0ZEd8Yw5gfioiNEfGdIfZvFRErIuKycntBRJy3pfNJ6i0Tu12AJG2hVcDMBu1Pd7qQNjsdeDYz72jVgJn5SkT8EjgJ+IcGXY4GpgLXltsXAjdFxKXlb5VK6mOGP0n9amNm/rpbk0fEtpm5rs1zfAD4OvCDNgx/LXBFROyTmUvq9p0EvAjcBZCZd0bE6xQ/Y3VpG2qR1EGe9pU0LkXEnuXp0RMj4vKIWBURL0bE+WWoGtz3gPLU5prydX1E7DJo/xHlWEdHxE0R8Sbwz+W+qRExLyLWRsTLEXF2RPw4IpaW+38vIt6OiFPq5oyIeDYifjLM1/hj4KPAz0b4rgeUp4aviYgJg+a9PCJeKee/OyIOHfSxnwHrKYLe4LEmAV8E5uX7f//zBuCvhqtDUn8w/EnqWxExsf7VoNvfA28Cfw78O/C98s+1MfYFfgVMoljZOhX4ODA/IqJurCuA31GEoyvKtquAPwG+BcwGPgd8ufaBzFwB3AjMqhvrCGBPYGCYr/hZ4KnMfH2oDhHxCeB2ih+OPyUzN0XENsCisq4zgT8FXgMW1UJtefr2ZurCH/AFYAfeO+VbczfwyYiYOky9kvqAp30l9audgHfqGyNir8xcOqjpfzOzdmPDrRExEzgB+M+y7fvAMuCYzNxQjvEw8ATweWDBoLGuz8zNNz5ExAEUQfDEzLy+bLsNeIEicNZcAdwSEXtn5jNl2yzggcx8ZJjv+Eng0aF2lit5N1OE2m8OWqk7GTgA+HhmLi77LgKeBL5DEQihCHjzImL/zHy8bDuJInA+WDfd74AAZgC3DlOzpB7nyp+kfrUKOKTB6+W6frfUbT8O7Dpo+yiKlbl3B60ePgsspQg6gy2o267tn19rKK8DXFTX7zbgOeAUgIjYgSKADrfqB7ALsHyIfZ+iCGFzM/P0ulO0RwEPAM/WrYjeUfed5lOE1JPKuranCLz1q34MqmOXBvsk9RFX/iT1q42ZeX8T/ervTt1AcYq3Zhpwdvmqt1vd9it127sAazLz7br21wZvZGZGxADw1YiYA5xIcfz96Qi1T6K4Lq+Rz5VjXN1g3zTgMBqsjAKbb+7IzLci4iaK8Pc94HhgW2Beg8/V6pjUYJ+kPmL4k1R1tWvy/rXBvvpVt6zbXgbsEBGT6gLg9AZjDVCcYj6S4rrC/8rMlU3UtuMQ+35IscJ3a0R8pu6O3RXA/cDfNvhcfZi8FviLiDiYIgQ+lJlPNPhcrY4VI9QsqccZ/iRV3W0U18c9UHfqtBm1lccvUl5DGBHbUtxosWZwx8x8ISJuAc4HPk3jZxTWexLYa4h971DcuPILihs5Pp2ZL5X7bqNYGXw+M18dYY6FFIHutPIz5w7Rb8/y/akm6pbUwwx/kvrVxIg4rEH7C4NCUDPmAL8BFkTElRSrfR+lCHBXZebtQ30wMx+NiPnAv5TX8S0DzgDeAt5t8JErgOspnqHXzE0TvwL+LCI+kJn/b7zMXBcRx1FcY7goIv4wM1+jOBV8GnB7RPwYeIbiBpk/AJZl5j8OGuOdiLgB+Ouy6bohaplBcZ3lY03ULamHecOHpH41Bbinwav+kSrDysynKK6PewuYC/wPxercepr7tZBTKcLXJcCVFDdV3AysbtD3v4GNwL81CnMN/JziGrxPDVP/m8AxZb0LI2JKeQr6SIqAeT7FTS//BOxHEXTrXUtxJ+/dmfn8EFPNBG5ssm5JPSxGf5ZDkjSU8s7aR4F7M7P+wc6fpwiAH8vMpn6GLiJ+DryYmV9vebFNiogpFDe7HJWZd3WrDkmtYfiTpDGIiC8BHwEeASYDX6NYJTs8M39T9vkIxarbpRTX4X1hFOMfQnEN3x5N3CDSFhFxDjAzM4/oxvySWsvTvpI0NmspTjXfRHH6dDpwXC34lWZTBLi3gdNHM3hm3gecBezekmq3zCrgm12cX1ILufInSZJUIa78SZIkVYjhT5IkqUIMf5IkSRVi+JMkSaoQw58kSVKFGP4kSZIq5P8AFLS0fkHhokwAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 720x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(10,8))\n",
    "plt.loglog()\n",
    "plt.errorbar(x081, y081, yerr=yerr081, xerr=xerr081, fmt='^', label='2008(1)')\n",
    "plt.errorbar(x101, y101, yerr=yerr101, xerr=xerr101, fmt='s', label='2010(1)')\n",
    "plt.errorbar(x102, y102, yerr=yerr102, xerr=xerr102, fmt='*', label='2010(2)')\n",
    "plt.xlabel(\"Energy (keV)\",fontsize=15)\n",
    "plt.ylabel(r\"$keV^2 (Photons{\\ }cm^{-2}{\\ } s^{-1}{\\ } keV^{-1})$\",fontsize=15)\n",
    "plt.legend(fontsize=13)\n",
    "plt.tick_params(labelsize=13)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
